A SPITZER/IRS SPECTRAL STUDY OF A SAMPLE OF 
GALACTIC CARBON-RICH PROTO-PLANETARY NEBULAE 

Yong Zhang and Sun Kwok 
Department of Physics, University of Hong Kong, Pokfulam Rd., Hong Kong 
zhangy96@hku.hk; sunkwok@hku . hk 
and 

Bruce J. Hrivnak 

Department of Physics & Astronomy, Valparaiso University, Valparaiso, IN 46383, U.S.A. 

bruce . hrivnak@valpo . edu 

ABSTRACT 

Recent infrared spectroscopic observations have shown that proto-planetary 
nebulae (PPNs) are sites of active synthesis of organic compounds in the late 
stages of stellar evolution. This paper presents a study of Spitzer/IRS spectra 
for a sample of carbon-rich PPNs, all except one of which show the unidenti- 
fied 21 /im emission feature. The strengths of the aromatic infrared band (AIB), 
21 /im, and 30 /im features are obtained by decomposition of the spectra. The 
observed variations in the strengths and peak wavelengths of the features sup- 
port the model that the newly synthesized organic compounds gradually change 
from aliphatic to aromatic characteristics as stars evolve from PPNs to planetary 
nebulae. 

Subject headings: infrared: stars — stars: AGB and post-AGB — stars: circum- 
stellar matter 



1. Introduction 

Proto-planetary nebulae (PPNs) represent a relatively short stage (~ 10 3 yr) when a 
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with B-G spectral type and are not yet hot enough to ionize the surrounding AGB remnant 
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envelope. As descendants of AGB stars, PPNs are expected to inherit many of the observed 
properties of AGB stars, including infrared excess and molecule-line emission from the dust 
and molecular components of the circumstellar envelopes of AGB stars. By comparing the 
spectroscopic properties of the dust component from AGB stars to PPNs to PNs, one sees 
evidence for a continuous evolution as the result of chemical synthesis and a hardening 
radiation field ( IKwokl 12004 ). While the infra red spectra of AGB stars are ch aracterized by 
silicates or silicon carbide emission features ( IKwok. Volk. fc Bidelmanlll997l ). the aromatic 
infrared bands (AIBs) dominate the infrared spectra of carbon-rich PNs. How this chemical 
transition occurs is one of the most interesting question in astrochemistry. PPNs, being 
in the evolutionary stage between these two phases, offer important clues to the chemical 
synthesis of organic compounds in the late stages of stellar evolution. 

Specifically, while carbon-rich PNs show strong emission features at 3.3, 6.2, 7.7-7.9, 
8.6, 11.3, and 12.7/zm due to aromatic compounds, these features are not seen in AGB 
stars and they only begin to emerge during the PPN phase. Furthermore , PPNs show 



aliphatic features at 3.4 and 6.9 /im which are generally weak or abs ent in PNs ( iGeballe et al. 



19921 ; iKwok. Volk. & Hrivnaklll99a Kwok. Volk. fc Bernathll200lh . Whether this aliphatic- 
aromatic transition is the result of chemical synthesis or radiative interactions is yet to 
be determined. I n a comprehensive study of the peak wavelengths of the AIB features, 
Sloan et al.l (120071 ) found that the peak wavelengths of the 7.7 and 11.3 /im AIB features are 
correlated with the temperatures of the central stars, which they interpret to be the result 
of change in the aliphatic to aromatic content ratio. 

The AIB features, although commonly assumed in astronomical literature to be due to 
polycyclic aromatic hydrocarbon (PAH) molecules, are in fact more likely (at least in the cir- 
cumstellar environment of late-type stars) to arise from complex, amorphous, impure organic 
solid-state compounds (IKwokl |2009| ). Examples of such amorphous compounds include hy- 

droge nated amorphous carb on (HAC;|Jones. E>uley. fc Williamslll990l ; IScott. Duley fc Jahani 

1997), carbon nanoparti cles jHu &: Dulevll2008l ). quenched carbonaceous composites (jSakata et al 
19831 ). coal or kerogen ( PapoularT2001 ). Common properties among these candidates are a 
mixed sp 2 /sp 3 chemistry, and in the case of coal or kerogen, the presence of elements such 
as O, S, N, in addition to C an d H. A model of such compounds has been presented by 
Pendleton fc Allamandolal (l2002f ). 



To complicate the issues further, it is now known that the infrared spectra of C- 
rich PPNs contain two unidentified broad emission features at around 21 and 30 /im (see 
Hrivnak et al.ll2009l . and references therein). The 21 /im feature was first discovered in spec- 
tra taken by the In frared Astronomical Satellite (IRAS) Low Resolution Spectrometer (LRS) 
(IKwok et al.lll989l ). Based on an Infrared Space Observatory (ISO) study of eight 21 /im 
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sources, IVolk et al.l (119991 ) found that they all have the same intrinsic profile and peak wave- 
length ( 20.1 /im). To date, the number of known 21 /im sources is very small, ~16 in the 



Gala xy (IHrivnak et al.ll2008l ) and ~8 in the Large and Small Magallenic Clouds (IVolk et al. 
2010l ). The carrier of the 21 /im feature remains unknown. A number of candidates, in- 



cluding TiC, AIBs, fullerenes, hydrogenated amorphous carbon (HAC), nanodiamonds etc, 
have been proposed (see Speck & Hofmeister (2004) and Posch et al. (2004) and references 
therein). The fact that all or almost all the 21 /im sources are C-rich PPNs suggests that the 
carrier of this feature is carbon-based and is associated with a very rapid chemical processing. 

The very broad 30 /im feature, usually extending from 25 to 45 /im, was first discovered 
by Forrest et al.l ( )198ll ) in C-rich AGB sta rs and PNs, and subsequently has been commonly 
detected in C-rich PPNs (jVolk et al.l 120021 ). An ISO study of a large sample of evolved stars 
shows la rge variations in t he peak wavelength and width of the 30 /i m feature among d ifferent 
sources Jlionv et alJbooj ). On the basis of the ISO spectra of PPNs, IVolk et al.l hoO$ \ found 
that the feature can be resolved i nto a broad component at 33 /im and a narrower component 
at 26 /im, although ISloan et al.l ( 120031 ) argued that the 26 /im feature is an artifact caused 
by an inaccurate calibration. The origin of th e 30 /im feature is still u n der debate. It was 
first suggested as solid magnesium sulfide (MgS; lGoebel fc Moseleylll985l ). iHony et al.l (120021 ) 
modeled the profile using MgS dust grains with a temperature different from the bulk of dust 
that contributes to the continuum. They argued that the 26 /im feature c an be explained 
by MgS with a different shape distribution. However, IZhang et al.l ( 120091 ) found that the 
available MgS dust mass in the circumstellar envelope is too low to account for the observed 
feature strength. Moreover, the fact that the feature is only detected in C-rich sources and 
takes up a significant fraction (~ 20%) of the total luminosity of the object seems to favor 
the more abundant carbonaceous material as possible carrier. 

The infrared spectra of C-rich PPNs therefore consist of contributions from several 
components: a base continuum commonly assumed to be due to thermal emission from 
amorphous carbon grains, a component due to aromatic and aliphatic compounds, and 
broad emission features at 21 and 30 /im. With the Spitzer spectroscopic survey of PPNs, 
we now have a large enough sample to derive the systematic behavior of the spectral features 
and draw conclusions on their possible evolution during the transition from the AGB to the 
PN phase. 

In this study, we present Spitzer spectra of a sample of carbon-rich PPNs, with the goal 
to investigate the spectral behavior of the AIB, 21 /im, and 30 /im features. In order to derive 
the intrinsic strengths and profiles of the broad and often overlapping emission features, we 
need to correctly subtract the continuum under the features. A procedure of optimized 
fitting is used to decompose the observed spectra into individual emission components. The 
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paper is organized as follows. In Section 2, we briefly describe the data used in the current 
study. The method for spectral decomposition is presented in Section 3. In Section 4, we 
report the measurement results and discuss the implications of our findings. The last section 
summarizes the conclusions. 



2. Data 

The sample includes ten PPNs, all of which except one (IRAS 01005+7910), are known 
21 /im sources. The observations were carried out with t he Infrared Spectro graph (IRS; 



Houck et al. 2004) on the Spitzer Space Telescope (Spitzer, IWerner et al.ll2004l ) at different 
epochs between 2004-2008. Most of the data were obtained from the Spitzer archive as part 
of the program No. 20208 (PI: B. Hrivnak) and 93 (PI: D. Cruikshank). These sources were 
observed with the Short-High (SH) and Long-High (LH) modules, covering the wavelength 
range from 9.9-19.6 /mi and 18.7-37.2 /im, respectively. The spectra of IRAS 07134+1005 
and IRAS 20000+3239 (with extremely strong 21 /zm emission) are not included in our anal- 
ysis because the 21-/im feature is saturated in the observations. We also include the data of 
three targets (IRAS 01005+7910, IRAS 04296+3429, and IRAS 22223+4327) observed with 
the short-Low (SL) module in Spitzer programs 30036 (PI: G. Fazio) and 45 (PI: T. Roellig), 
which cover a larger wavelength range (5.5-37.2 /mi). These allow the study of the shorter 
wavelength 6-10 /xm AIB features. In some cases there are discontinuities in the fluxes be- 
tween the different modules. This is due to different slit sizes and slight mis-pointing. To 
ensure that all spectra were on the same flux scale, the regions of spectral overlap are com- 
pared and used to scale and combine the spectra derived with different modules for each 
object. In the case of IRAS 19441+2401, the star is outside of the SH aperture so the 
spectrum would only have recorded part of the extended emission. Consequently, the short- 
wavelength part of t he spectrum m a y not be reliable. Details of the data processing are the 



same as described in iHrivnak et al.l (120091 ) . The extracted spectra are displayed in Figure [TJ 



3. Model fitting 

The IRS spectra of all ten sources show strong infrared excess due to dust continuum 
emission. Superimposed on the continuum are several strong, broad emission features. There 
is no evidence of strong atomic emission lines as one would expect in a PN. The only narrow 
features present are molecular lines. This is consistent with the fact that the central stars of 
these objects are not hot enough to photoionize the surrounding circumstellar materials. In 
this section, we report on our efforts to fit the observed spectra with a combination of dust 
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continuum, dust features, and H 2 lines. 

For the spectral d ecomposition, we used the IDL package PAHFIT originally developed 



by ISmith et al.l (120071 ) to fit the Spitzer IRS spectra of nearby galaxies. The model spectra 
take into account the contributions from stellar continuum, thermal dust continuum, H2 
lines, and AIBs. We modified this code to include fittings to the 21 and 30 /im features. 
The optimal fitting to the observed spectra is achieved through the Levenberg-Marquardt 
least-squares algorithm. The fitting results are shown with the spectra in Figure [TJ 

In our model, the stellar continuum is assumed to be a 5000 K blackbody for all the 
sources. While the individual stellar temperatures are known to be mostly in the range 
of 5250-7250 K, the actual stellar contribution to the mid-infrared spectra is minimal (see 
Fig.[T]). In almost all cases, the effects of photospheric continuum contributions are negligible. 

The dominant contribution to the mid-IR continuum is due to dust thermal emission. 
A modified blackbody model [I\ ~ \~ 2 B\(T), where B\(T) is the blackbody function with 
a temperature T] is adopted to simulate the dust thermal continuum. However, for most 
of the sources, we must use two components with different temperatures to get the best fit. 
The warm dust component has temperatures from 140 to 180 K, and the cold component has 
temperatures from 60 to 80 K. The only exception is IRAS 19477+2401, which does not need 
a warm component (although we note that the short-wavelength spectrum of this object may 
not be reliable). 

Eight H2 rotational lines from S(0)-S(7) are included in the fitting. The line profiles are 
assumed to be Gaussian. The H2 0-0 S(0) para ground-state transition at 28 /im is clearly 
detected in IRAS 04296+3429 and IRAS 22223+4327. From the spectra, we can see that the 
H 2 lines are much narrower than dust features. 



Also included in the fitting are all the AIB features in ISmith et al.l (120071 see their 

2 

Table 3). Drude profiles [I\ ~ Q^p^z~\^jj^FipyZ 1 where Ao is the central wavelength, and 7 is 
the fractional FWHM] are assumed for each feature. Figure [1] shows that the 7.7, 11.3, and 
12.7/xm complexes can be well fitted with blended subfeatures. One should bear in mind 
that these subfeatures do not necessarily represent distinct physical emission bands. 



The intrinsic shape of the 21 /im feature has been investigated by IVolk et al.l ( 119991 ). 
who found that the profiles are similar among different sources, all showing "a sharp rise 
on the short-wavelength side, beginning at 18 .5 /im to a peak at 2 0.1 /im, and a gradual 



decline on the long-wavelength side to 24 /im" . iHrivnak et al.l (120091 ) found the same shape 
in the Spitzer spectra they studied and we see similar shapes for the three additional objects 
(IRAS 01005+7910, 04296+4329, 22223+4327) in our sample (see Figured). To represent 
this profile asymmetry, we have assumed an asymmetrical Drude profile with the 7 value in 
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the long-wavelength side three times larger than that in the short-wavelength side. 



The very broad 30 /im feature is obviously present in the Spitzer IRS spectra. iHony et al. 



( 120021 ) found that the peak wavelength of this feature shifts from 26 /iin in some AGB stars 
to 38 /iin in PNs. Through the fitting, we find it difficult to reproduce the observations using 
a single component because the decomposed profile of this feature differ from star to star. 
Although the Spitzer spectra do not show the two distinct components in the 30 fim feature 
previously deduced from the ISO spectra, we find it practical to use two Drude profiles 
peaking at 26 /im and 30 /iin to represent this feature as they provide better fits to the data. 

Dust extinction is not considered in the present study as its effect is expected to be small 
in mid-infrared wavelengths. If present, extinction may reduce the emission flux at 9.7 /im 
and 18 /im due to interstellar silicate absorption. In some cases, our model cannot reproduce 
the fluxes around 9.7 /im and 18 /im (the model fluxes at these wavelengths are slightly higher 
than observed), which might be an indication of the presence of interstellar dust extinction. 
Nevertheless, the small discrepancy does not affect our study of dust features. 

Although this study is mainly focused on individual features, we explore the possibility 
of using principal component analysis (PCA) to study global spectral properties and classify 
PPNs. Our results are presented in Appendix A. 



4. Results and discussion 

We are able to obtain reasonably good fits to the observed spectra using the above 
model and the results are shown in Figure [TJ In Table [T] we give the integrated fluxes of 
the continuum and the 21 and 30 /im features. The flux fractions of the 21 and 30 /im 
features to the total IR spectra are also presented in this table. Table [2] lists the strengths 
of all the major AIB features. Since the SL observations were made only for three objects 
(IRAS 01005+7910, IRAS 04296+3429, and IRAS 22223+4327), we are only able to obtain 
the 6.2, 6.9, 7.4, 7.8, 8.3, and 8.6 /im fluxes for these three objects but not the others. 
The feature strengths are derived by integrating the fluxes above the model (star + dust) 



conti nuum. To extract the strengths of dust features, some authors (e.g.. lBernard-Salas et al. 



20091 ) derived the continuum through fitting a spline function to selected spectral regions 
free from feature emission. From the PPN spectra shown in Figure [TJ we can see that the 
features are broad and are blended with each other. Consequently, it is difficult to find a 
region completely free of dust features to adequately define the baseline. Given the difficulty 
in selecting appropriate sections of the spectra to represent the continuum, we believe that 
our present method gives more robust results. 
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For almost all of the targets, two dust components with different temperatures are re- 
quired to reproduce the thermal continuum. The blackbody temperatures used for the fitting 
are listed in Table [TJ The cooler component probably arises from an extended circumstellar 
envelope left over from the AGB progenitors, and the warm component may originate from a 
central torus or more recent mass loss. Some PPNs have been found to exhibit equatorial tori. 
For example, the dust con tinuum and CO map of the 2 1 /im source IRAS 07134+1005 sug- 



gests the presence of torus ( Kwok. Volk. fc Hrivnakll2002l ; lMeixner et al.ll2004j ; iNakashima et al. 



2009|). 



Tables [T] and Table |2] also give the fitting uncertainties of the feature strengths, which 
are estimated using the full covariance matrix of the least-squares parameters. We can 
see that the errors are typically less than 20%. Adding more blackbody components with 
temperature of 50-200 K to the fitting might lead to 5 — 15% changes in the fluxes of the 
broad 30 /im feature, but will have minimal effect on those of the other features (< 5%). 



4.1. 21 /im and 30 |im features 

The strengths of the 21 /im feature vary among the sources. In the strongest 21 /im 
source (IRAS 06530— 0213) in our sample, the fraction of the flux emitted by the 21 /im 
feature to the total emitted infrared flux amounts to ~5%. Since we fixed the peak wave- 
length and shape of the 21 /im feature in the fitting process, the good agreement between the 
model and observation supports the su ggestion that the intrinsic profile of the 21-/xm feature 



does not change in different sources (jVolk et al.l If 999l ). We also note that the two PPNs 



(IRAS 19477+2401 and IRAS 22223+4327) with weak 21 /im features have no or weak warm 
dust component. This might indicate that the 21 /im feature mostly arises from the warmer 
torus, although infrared imaging is needed to confirm this suggestion. The one object in 
our sample without the 21 /im feature is IRAS 01005+7910. It has a stellar temperature of 
21 000 K, which is considerably hotter than the other PPNs (see Tabled]). The presence of 
the 21 /im feature in PPNs but not in AGB stars or PNs has previously been interpreted as 
an indication of the evolution of the feature due to chemical synthesis or changing excitation 
conditions. 

The 30 /im feature is present in all of the PPNs. Its flux typically accounts for ~20% 
of the total infrared flux in the PPNs (they range from 8 to 28%). From Figured] one can 
see that the 30 /im feature has an asymmetric profile with a sharp blue wing and a broad 
red tail. Although the presence of two se parate components for this feature as seen in the 



ISO spectra has been shown to be unreal (IHrivnak et al.ll2009l ). we find it helpful to use two 



components to fit this broad feature because a single component of fixed shape and central 
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wavelength cannot give a satisfactory fitting. We should note that the two components do 
not necessarily mean two individual features. We find that the relative flux ratios of the 
two components vary among so urces. This r epres entation of the changing shape and peak is 
consistent with the findings of iHony et al.l (120021 ) who simulated the MgS emission feature 
and found that the peak position changes with the grain shape and temperature. 



Bernard-Salas et al.l (120091 ) investigated the Spitzer/IRS spectra of 25 PNs in the Mag- 



ellanic Clouds (MCs) and found that about half of them possess the 30 /im feature. They did 
not find a clear relation between the presence of this feature and nebular sulfur abundance, 
which is used to argue against the suggestion of MgS as possible carr i er. T he strengths of 
the 30 /im features in MC PNs (see Figure 6 in iBernard-Salas et al.l 120091 ) are similar to 
those in MW PNs, suggesting that the different metallicity of the two galactic systems is not 
a major factor in the existence of the 30 /im feature. However, when comparison is made to 
the strengths of the feature in MW PPNs as presented in this study, one sees that the feature 
is stronger in PPNs. This suggests that the 30 /im carrier is gradually being destroyed or 
not excited in the PN phase. If MgS is the carrier, we expect the amount of sulfur released 
into the gaseous nebula will in crease as the star evolve to a PN. However, the low sulfur 



abundance found in PNs (e.g. iPottasch fc Bernard-Salad 120061 ; IBernard-Salas et al.l 120081 ) 



does not seem to support this idea (although there exist other explanations for the 'sulfur 
anomaly'; see Henry et al. 2006). The decline of the strength of the 30 /im feature is more 
likely to be the result of changing physical or chemical conditions. 



4.2. AIBs 

The AIB features are prominently present in all of the PPNs in this sample and the peak 
wavelengths of the features detected are given in Table [3j The profiles and central wave- 



lengt hs of AIBs are known to vary in different astronomical sources (see, e.g. iPeeters et al. 



20021 ). Almost all of the sources have the 11.3 /im feature peaking at the same wavelength, 
11.39+ 0.03 /im. The exceptions are IRAS 22574+6606 and IRAS 01005+7910, the latter 
of which is a much hotter source and does not possess the 21 /im feature. 

In order to examine the features in detail, we have plotted in Figure [2] the continuum- 
subtracted spectra of the three PPNs observed in the IRS mode SL 6-10 /im region together 
with the spectra of the young PN IRAS 21282+5050. The peak wavelengths of the aromatic 
and aliphatic emission features are marked. The 6.2 /im aromatic feature is detected in all 
three PPNs that were observed in this wavelength region. The 6.9 and 7.25 /im aliphatic 
features are particularly distinct in IRAS 04296+3429, whereas the 8.6 /im aromatic feature 
is relatively weak in that object. In IRAS 22223+4327, the 6.9 /im is not obvious from 
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looking at the spectrum alone, but from the fitting results, the integrated strength of this 
feature is quite high due to its width. 

One can com pare the 6.2 /im/11.3 pm in tensity ratio to that found in PNs in the MW, 
LMC, and SMC (IBernard-Salas et alJ EooiL see their Fig. 2). For IRAS 04296+3429 and 
IRAS 22223+4327, the ratio is low, 0.45 and 0.47, respectively, while for the hotter PPN 
IRAS 01005+7910, the ratio of 1.40 is similar to that found for most of the PNs. An 
examination of the shorter wavelength features also shows shifts in the peaks of the features 
between IRAS 01005+7910 and the other two sources observed and also the appearance 
of a small shift in the peak of the broad 7.7-8.2 /im complex between IRAS 22223+4327 
and 04296+3429. Also seen is a difference in the peak at 7.60 /im for the hot PPN IRAS 
01005+7910 and 7.85 /im for the PN IRAS 21282+5050. There seems to be a general trend 
for the peak wavelength of the 7.7 /im feature to shift to shorter wavelengths with evolution. 

In Figure [3j we have plotted the peak wavelengths of the 11.3 /xm feature as functions 
of the fractional strengths of the of the 21 /im feature. The only two objects with redshifted 
11.3 /xm are the ones with weak or no 21 /im feature. The peak wavelengths of the 11.3 /im 
feature as a function of central star temperatures are plotted in Figure |H For co mparison, we 
have also included in Fig ure H] the values for Ae/Be and other C-rich stars from lSloan et al. 



(120071 ) and iKeller et al.l (120081 ). Although these objects do not belong to a homogeneous 
sample, this figure seems to suggest that the peak of the 11.3/im feature occurs at shorter 
wavelengths at higher central-star temperatures. It has been suggested that the sources 
with most redshifted AIB fe atures are exposed to weaker UV radiation fields and are less 
processed (ISloan et al.l 120071 ). While we note that the amount of UV radiation background 
in our PPN samp le is minimal, our results are consistent with the trends as reported by 
Sloan et al.l (120071 ). Proposed explanations for the wavelength change of the AIB features 
include nitrogen-inclusion, variations in the fraction of carbon isotope, size distribution of 
the grains, ionization stat e, and the content ratio of aliphatic to aromatic bands (see, e.g., 
Bernard-Salas et al.ll2009l . and references therein). We suggest that this shift of the peaks of 
the AIB features to shorter wavelengths is related to the relative strengths of the aliphatic to 
aromatic compone nts as well as growth of the siz es of the aromatic rings as the star evolve 
from PPN to PN (IKwok. Volk. fc Bernathl l200ll ). Photochemistry may be the underlying 
cause of these changes. 

The 11.3 fim AIB feature is identified with the C— H out-of-plane bending mode of 
aromatic compounds. However, when the number of rings in the aromatic units is small, 
there can be 2 (duo), 3 (trio), or 4 (quarto) adja cent CH groups and their out -of-plane 
bending mode frequencies can be slightly different (jHudgins fc Allamandolalll999l ). Gener- 
ally, the duo-CH bending mode is in the 11.6-12.5 /im range, the trio-CH mode is in the 
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12.4-13.3 /im, and the quarto-CH mode in the 13-13.6 /mi range. While the 11.3 /im solo- 
CH mode is dominant in PNs , the other modes are seen to be relatively stronger in PPNs 
(IKwok. Volk. fc Hrivnaklll999l ). This observation has l ed to the sug gestion that the aromatic 
units grow in size as the star evolve from PPNs to PNs (lKwokll20041 ). With our present fitting 
exercise, it is possible to quantify this change. In Figure |5] we plotted the ratios of the sum of 
the strengths of the duo, tri, and quatro modes to the solo mode as a function of central star 
temperature. It can be seen that the strengths of the duo, trio and quatro modes decrease 
relative to the solo mode as the star evolves to higher temperatures, supporting the scenario 
that the aromatic units get larger with evolution. 

The derived fluxes from each of the AIB features as well as their total contribution to 
the observed infrared fluxes are given in Table [2J Comparing the last column of Table [2] 
to the values in columns 8, 10, 12 in Table (U we can see that the contribution of the AIB 
features to the total emitted infrared flux is much smaller than that of the 30 /im feature 
and is of the same order (a few %) as that of the 21 fim feature. In some cases, the flux of 
the 11.3 /jm feature alone is larger than (e.g., in IRAS 05341+0852 and IRAS 07430+1115), 
or about equal to (as in IRAS 22574+5435) that of the 21 /im feature . 



Cerrigone et al.l (120091 ) found that about 40% of their post-AGB star sample exhibit 



both AIBs and silicate features, and suggested that the spatial segregation of different dust 
populations in the envelopes can lead to the mixed chemistry and a disk/torus can serve as 
a reservoir preserving O-bearing molecules. However, we do not detect silicates features in 
our (C-rich) sources, suggesting that the 21 /im sources do not favor the format ion/ survival 
of O-bearing compounds. 



4.3. 15-20 /im plateau feature 



A 15-20 /im plateau emission feature is commonly seen in PNs and other objects with 
a hot central star and has been attributed to C-C-C bending modes of aromatic rings 
( IVan Kerckhoven et al.l 120001 ). For example, this fea ture is seen in IRAS 06556+1623 and 
IRAS 18442—1144 and several other post-AGB stars ( jCerrigone et al.ll2009l ). In most of our 
sources, this region is difficult to examine because of overlap with the blue wing of the 21 
/im feature. However, in IRAS 01005+7910 this is seen as a prominent feature, as shown in 
Figure El Integrating the fluxes between 15.2 and 20 /im, we obtain an integrated strength 
of 8.5 xlO -13 W/m 2 , leading to strength ratios relative to the C-C stretching mode at 6.2 /im 
of 22 (Ipiateau/ h.2) and the C-H out-of-plane bending model at 11.3 /im of 31 (I p i a teau/ hi.3)- 
The plateau/6.2,11.3 strength ratios are much large r than the values (<5) found i n the sam- 
ple of H II regions, Ae/Be stars, and PNs studied by lVan Kerckhoven et al.l (120001 ). Whether 
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this is an evolutionary trend or particular to IRAS 01005+7910 is difficult to say based on 
only one case. In this respect, it would be good if one could obtain 15-20 /on spectra of 
other C-rich PPNs without a strong 21 /im feature. 



4.4. 15.8 /im feature 



In a previous paper (IHrivnak et al.l 120091 ). we called attention to a moderately broad 
emission feature at 15.8 /im seen in the Spitzer spectra of several carbon-rich PPNs with 
the 21 /im feature. This feature is seen in the several additional spectra included in this 
study. While a weak feature at 15.9 /im has been detected in the refl ection nebula NGC 7023 
( iWerner et al.l 120041 ) and some star-forming galaxies (see Fig. 5 of ISmith et al.l 120071 ). the 
feature seen in the PPNs in the present study is significantly broader. We noted there 
was a suggestion of a correlation with the strength of the 21 /im feature in the sense that 
the two sources with the strongest 21 /im features also had the strongest 15.8 /im features 
(IRAS 06530—0230 and 23304+6147). We can examine this more quantitatively in this study 
with the measurement of the strength of this feature in these sources (see Table 2) and the 
inclusion of two other 21 /mi sources (IRAS 04296+4329 and IRAS 22223+4327). 

Fig. [7] displays the 15.8 /im fluxes versus the 21 /im fluxes for the PPN sample. With 
the exception of IRAS 04296+3429, there seems to be a positive correlation between the two 
fluxes, although more data, specially those of strong 21 /im sources, are required to confirm 
this conjecture. Both 15.8 and 21 /im features have an asymmetrical profile. In Fig. |8l 
we compare the profiles of the two features. While in all the cases the 15.8 /im feature is 
narrower than the 21 /im feature, they are similar in shape: both have a steeper blue rise 
and a shallower red tail. 



4.5. Molecular lines 

In addition to t he ab ove broad features, there are also molecular lines in the spec- 



tra. 



Hrivnak et al.l ( 120091 ) have detected the 13.7 /im z/ 5 line of C2H2 in five PPNs, one 
(IRAS 22574+6609) in absorption, three (IRAS 05341+0852, IRAS 06530-0213, and IRAS 
07430+1115) in emission, and one (IRAS 23304+6147) in both absorption and emission (P 
Cygni profile). From our examination of the spectra, the C2H2 could also be present in 
emission in IRAS 04296+3429. C2H2 is the most abundant molecular species after H2 and 
CO in carbon-rich envelopes. In our sample, about 70% PPNs exhibit the C2H2 feature. 
IRAS 05341+0652 has the strongest C2H2 emission. We do not find correlation between the 
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fluxes of C2H2 and those of other features. 

Furthermore, we can unambiguously identify the H 2 0-0 S(0) line at 28 /im in IRAS 04296+3429 
and IRAS 22223+4327 and their spectra are shown in Figure. M However, the expected next 
strongest H 2 line ortho ground-state v = — S(l) line at 17 /im is not detected. We do not 
find any correlation between the strengths of dust features and those of molecular hydrogen 
emission. IRAS 04296+3429 shows relatively strong AIB and 21 fim features, but these fea- 
tures are weak in the spectrum of IRAS 22223+4327. Table @] gives the fluxes of the S(0) 
transition and the flux upper limits of the S(l) transition estimated from the 3cr flux noise 
level. These values allow us to determine the upper limits of excitation temperatures in the 
two objects. Under an assumption of local thermodynamic equilibrium (LTE), we have 



Ny,j = N exp(-^j/fcT cx ) 
9 j tot S 1 ,/ iJ /5(j/exp(-E, ; / i j//A;T ex )' 
where N v> j is the population of the upper level (v, J), gj is the level degeneracy under the 
assumption of ortho/para ratio = 3, (gj = (2 J + l) for even J and 3(2J+ 1) for odd J), N tot 
is the total H 2 density, E Vj j is the energy of the upper level, T ex is the excitation temperature, 
and k is Boltzmann's constant. The flux of a given transition, F, is proportional to N v j 
(F ex N Vt jA/X, where A is the radiative transition probability). From Equation (JT|) we have 



E v J Jk - E v j jk 



ex ~ HFoXog^A,) - HF^g^Ao) ' 1 ] 

where the subscripts, and 1, represent the t ransit ions S(0) and S(l), respectively. Taking 
the E v ^j and A values from iRosenthal et al.l (120001 ). we derive the T ex upper limits of 56 K 
and 51 K for IRAS 04296+3429 and IRAS 22223+4327, respectively. The extremely low 
excitation temperatures of H2 seem to suggest that the contributions from UV pumping and 
shock waves to the excitation of H 2 are insignificant in the two PPNs. 

This is consistent with previous observations o f the H 2 v = 1 — S(l) line at 2.12 /im of 
these two and other PPNs ( Kelly fc Hrivnak 20051 ). It was found that H 2 emission was seen 
in F-G spectral type PPNs only if they were clearly bipolar in morphology with a relatively 
dense equatorial dust torus and a bipolar outflow, and that the H 2 emission in these cases 
was collisionally excited. 



5. Concluding remarks 



We have investigated the mid-IR spectral features - AIBs, 15-20 /xm plateau, 15.8 /im, 
21 fim, and 30 /mi - in a sample of PPNs using their Spitzer/IRS spectra. The strengths of 
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these emission features are determined using a decomposition method. We find that the 21 
and 30 /iin features are strong in these PPNs, contributing up to 5 and 30%, respectively, of 
the total infrared fluxes observed. A prominent 15-20 /mi broad emission plateaus is detected 
in IRAS 01005+7910. The 15.8 /mi emission feature, seen in almost all our objects, is found 
to have similar intrinsic emission profile as the 21 /mi feature. Molecular hydrogen lines are 
observed in two of the PPNs, showing the presence of the neutral circumstellar component. 
Most of the PPNs require two dust components to fit their underlying continua. We suggest 
that the cold component is due to an extended dust halo probably left over from the AGB 
phase, and a warmer component is due to a dust torus in the central region of the PPN. 
This suggestion can be tested by future infrared imaging observations. 

The fact that the duo, trio, and quarto C-H out-of-plane bending modes decline in 
strength with evolution relative to the solo mode is consistent with the suggestion that the 
aromatic component grows in size as the star evolves to higher temperatures. There is also 
evidence for the gradual shift of the peak of the 11.3 /im feature to the blue with evolution. 
The 21 and 30 /im features, which carriers are unknown, are shown to decline in strength 
with evolution. These observations will help us to constrain future chemical evolution of 
organic synthesis in the late stages of stellar evolution. 

This work was supported by a grant to SK from the Research Grants Council of the 
Hong Kong Special Administrative Region, China (Project No. HKU 7020/08P) and a grant 
to YZ from the Seed Funding Programme for Basic Research in HKU (200909159007). 



A. Principal component analysis 



We also investigated the variety of the PPN spectra using the method of principal 
component analysis (PCA). PCA is a mathemat i cal tool which has been commonly applied 



to spectral classification (e.g. iFrancis et al.lll992l ; ISuzukill2006l and references therein). The 



basic idea of PCA is to obtain new independent variables, called the principal components 
(PC), through orthogonal transformation. The new variables are a linear combination of 
the original data. The relative contributions to the spectrum variance of each PC can be 
determined. The PCA technique can reduce the multidimensionality of t he data and enables 



one to study the correlations between different spectral properties (see, IFrancis et al.lll992 



for a more detailed description of PCA). This analysis has the potential to give new insight 
into the spectral features. 



In our work, the common wavelength coverage (10-35 /mi) of the ten PP N spectra was 
considered for PCA. The formulation of PCA can be found in many papers (e.g. iRonen et al. 
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19991 ). We first subtracted the continuum and normalized each spectrum by its integrated 



flux. The mean spectrum is obtained by 



1 - 

<F(A,)>=-^F,(A,), (Al) 

i=l 

where N = 10 is the number of the PPN spectra, and j = 1, ...,M is the number of the 
wavelength channels. The data set can be denoted by a, N x M matrix 

X ij = F i (\ j )-<F(\ j )>. (A2) 

The eigenvector and eigenvalue of its covariance matrix, respectively, give the new axis 
(PCs) and the variance along the new axis. The PCs are ordered according to their relative 
contributions to the spectrum variance, such that the first PC accounts for the maximum 
variance. 

Figure [10] shows the mean spectrum and the first three PCs (PCI, PC2, and PC3). 
The first three PCs account for 79% of the variance in the spectra. The mean spectrum 
mainly exhibits three features: the AIBs around 10-15 /zm, the 21 /xm feature, and the 30 /im 
feature. PCI accounts for 37% of the variance. Compared to other PCs, the most prominent 
feature of PCI is the dip around 18 /im. Thus, it can modulates the strengths of the 21 /im 
feature and the emission from 12-15 /im. PC2 accounts for 29% of the variance. It displays 
the modulation of the AIB features from 10-15 /xm, and especially the 11.3 /im feature. 
Since PC2 can also modulate the flux of 12-15 /im emission, PCI particularly displays the 
modulation in the flux of the 21 /im feature. PC3 account for 13 percent of the variance. It 
reflects the variation in the relative fluxes of individual AIB features. We can conclude that 
the variance in this spectrum sample is mainly caused by the 21 /im and AIB features, but 
not by the 30 /im feature. 

Figure [TT1 shows the projections of the PPN spectra on to the first two PCs (denoted by 
pel and pc2). The values of pel and pc2 represent the contribution of each component to 
the difference between individual PPN spectrum and the mean spectrum. We can find that 
the non-21 /im source, IRAS 01005+7910, has the lowest pel value, and the strongest 21 /im 
source (IRAS 06530-0213, IRAS 23304+6147, IRAS 04296+3429) have the largest pel values. 
If the 21 /im feature arises from the warm torus(see Section 4), the value of pel could repre- 
sent the content of warm dust in the torus. The value of pc2 reflects the strengths of AIBs 
with respect to the mean spectrum. The sources with weak AIB features (IRAS 22223+4327 
and IRAS 19477+2401) have the largest pc2 values while the sources with strong AIB fea- 
tures (IRAS 05341-0852, IRAS 22574+6609) have the lowest pc2 values. Therefore, PCA 
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allows us to objectively and quantitatively classify the spectra of PPNs. However, the cur- 
rent sample is too small for the analysis of higher-order PCs. With a larger PPN sample, 
it would be intriguing to use PCA to study the global spectral behavior and the relation 
between the fluxes of different AIB features. 
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Fig. 1. — Decomposition of the IRS spectra. The black squares are the observations. The 
red solid lines show the thermal dust components. The stellar continuum, although included 
in the plots, are too weak to be seen. The gray lines show the total of the stellar and dust 
continuum components. The dust features and H 2 lines are indicated by the blue and violet 
lines, respectively. The green solid lines are total model spectra. 



- 19 - 



1 9477+2401 
30 I 




Wovelength (/im) 



22574+6609 

50 h 




Wavelength Qum) 



22223+4327 

80 




Wavelength (/J-m) 



23304+61 47 
100 




Wavelength (/xm) 



Fig. 1. — continued. 
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Fig. 2. — Continuum-subtracted IRS spectra. The continuum is the sum of the model stellar 
and dust continuum components. The solid and dash lines indicate the positions of aromatic 
and aliphatic emission features, respectively. 
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Fig. 3. — Peak wavelengths of the 11.3 jum feature versus flux fractions of the 21 /xm feature. 
The object with zero 21 jum flux is IRAS 01005+7910. 
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Fig. 4. — Peak wavelength of the 11.3 (ttrn f eature versu s effect ive temperatures of the central 
stars. The filled squares are objects from ISloan et al.l ( 120071 ) . the open squares are Ae/Be 
stars from lKeller et all ( 2008 ) and the filled circles are PPNs in the present study. 
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Fig. 5. — Ratio of the sum of the strengths of the 12.0, 12.6, and 13.5 /im features to the 
strength of the 11.3 /xm feature as a function of central star temperature. The filled circles 
are PPNs in our sample, and the filled squares are (from right to left) PNs IRAS 21282+5050, 
BD+303639, and NGC 7027. 
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Fig. 6. — The 15-20 /im emission plateau. The dashed line indicates the continuum. The 
continuum-subtracted spectrum is shown in the lower right panel. 
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Fig. 7. — The strength of the 15.8 /xm feature versus the strength of the 21 fxm feature. 
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Fig. 8. — The profiles of the 15.8 and 21 jum features in three PPNs. 
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Fig. 9.— The observed profiles of the H 2 0-0 S(0) transitions in IRAS 04296+3429 and 
IRAS 22223+4327. The dashed lines represent a Gaussian fit. 
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Fig. 10. — The mean spectrum and the first three PCs. 
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Fig. 11. — Projections on to the first and the second eigenvectors of the spectrum sample. 



Table 1. Effective temperatures and integrated intensities of the dust continuum and the 21 and 30yum features. 
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a Bascd on published model atmosphere studies. 
b The numbers in brackets are fitting uncertainties. 

c The flux fractions of the features to the total IR emission, which includes the contribution from the dust features from 10-35 ^tm and the dust thermal continuum. 
d G0 spectral type. 
e Al spectral type. 
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8.2 


)4296 - 


1-3429 


3.71(0.12) 


4.82(0.19) 


17.70(3.00) 


8.70(0.19) 


6.68(0.18) 


2.29E(0.14) 


1.09(0.23) 


8.26(0.27) 


8.54(0.34) 


6.21(0.31) 


3.15(0.25) 


0.27(0.02) 


9.1 


)5113- 


1- 1347 


b 












0.18(0.02) 


2.25(0.28) 


3.28(0.33) 


0.58(0.02) 


0.00(0.00) 


0.28(0.02) 


> 3.1 


)5341 - 


f 0852 














0.23(0.02) 


3.62(0.28) 


3.00(0.34) 


2.23(0.33) 


1.47(0.27) 


0.30(0.02) 


> 5.9 


)6530 - 


- 0213 














0.00(0.00) 


3.57(0.26) 


4.13(0.34) 


3.84(0.31) 


2.29(0.26) 


1.18(0.16) 


> 4.3 


)7430 - 


1- 1115 














0.99(0.02) 


4.83(0.28) 


5.79(0.34) 


2.71(0.33) 


1.42(0.26) 


0.39(0.02) 


> 4.2 


9477 -+ 


- 2401 d 














0.15(0.03) 


0.57(0.04) 


0.53(0.04) 


0.69(0.03) 


0.26(0.03) 


0.00(0.00) 


> 1.0 


>2223 - 


f 4327 


0.29(0.01) 


0.54(0.02) 


1.79(0.42) 


0.63(0.02) 


0.55(0.02) 


0.17(0.01) 


0.00(0.00) 


0.63(0.04) 


1.08(0.33) 


0.70(0.03) 


0.75(0.03) 


0.21(0.01) 


1.3 


!2574 - 


1-6609 














0.00(0.00) 


5.55(0.26) 


2.71(0.34) 


3.68(0.33) 


1.10(0.27) 


0.19(0.02) 


> 2.2 


!3304 - 


f6147 














0.00(0.00) 


2.64(0.15) 


9.49(0.32) 


6.87(0.32) 


4.31(0.26) 


2.12(0.14) 


> 2.6 



'In unit of 10~ 14 W/m 2 . The numbers in brackets are fitting uncertainties, 
'no feature strengths derived for the 7 objects without SL observations. 

: The flux fractions are the ratios of the sum of all the AIB features to the total IR emission, which includes the contribution from the dust features from 10-35 lira and the dust thermal 
itinuum. 

'The flux measurements of this object may not be reliable (see section 2). 



-32 - 



Table 3. Peak wavelengths of AIB features." 



IRAS ID 


6.2 (im 


7.7-8.2 /im 


8.6 /im 


11.3 /im 


01005 


+ 7910 


6.20 + 0.03 


7.57 + 0.05 


8.56 + 0.04 


11.25 + 0.02 


04296 


+ 3429 


6.26 ± 0.02 


8.02 ± 0.08 




11.42 + 0.09 


05113 


+ 1347 








11.41 + 0.13 


05341 + 0852 








11.40 + 0.09 


06530 


- 0213 








11.38 + 0.02 


07430 


+ 1115 








11.40 + 0.13 


19477 


+ 2401 








11.38 + 0.13 


22223 


+ 4327 


6.29 ± 0.03 


7.87 + 0.04 


8.53 + 0.03 


> 11.30 


22574 


+ 6609 








11.31 + 0.04 


23304 


+ 6147 








11.40 + 0.06 



a In unit of /im. 



Table 4. H 2 lines in IRAS 04296+3429 and IRAS 22223+4327. 



IRAS ID 


F(0-0 S(0)) 


F(0-0 S 


(1)) 






(W/m 2 ) 


(W/m 


2 ) 


(K) 


04296 + 3429 


3.2E-14 


< 4E- 


•16 


< 56 


22223 + 4327 


3.9E-14 


< 2E- 


-16 


< 51 



